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ABSTRACT 

Context. High- velocity stars in the Galactic halo, e.g. the so-called hyper- velocity stars (HVS), are important tracers 
of the properties of the dark matter halo, in particular its mass. 

Aims. A search for the fastest stars among hot subdwarfs (sdB) in the halo is carried out to identify HVS, unbound to 
the Galaxy, and bound population H stars in order to derive a lower limit to the halo mass. 

Methods. Based on the SDSS DR6 spectral database we selected stars with high rest-frame velocities. These radial 
velocity measurements were verified at several telescopes to exclude radial velocity variable stars. Out of 88 stars 
observed in the follow-up campaign 39 stars were found to have constant radial velocities. For twelve of them we 
measured a proper motion significantly different from zero and obtained spectroscopic distances from quantitative 
spectral analysis to construct the full 6D phase space information for a kinematical study. 

Results. All but one programme sdBs show halo characteristics, but can be distinguished into two kinematical groups, 
one (Gl) with low Galactic rotation typical of halo stars and a second one (G2) with rapid retrograde motion. We also 
investigate the possibility that the programme stars are not genuine halo stars but ejected from the Galactic disc or 
bulge. The Gl objects crossed the Galactic plane in the central bulge, whereas the G2 stars did in the outer Galactic 
disc. J1211-f 1437 (G2) is a HVS candidate, as it is unbound to the Galaxy if the standard Galactic potential is adopted. 
Conclusions. We conclude that in the ejection scenario Gl stars might have been formed via the slingshot mechanism 
that invokes acceleration by tidal interaction of a binary with the central supermassive black hole. The G2 stars, 
however, would originate in the outskirts of the Galactic disc and not in the central bulge. J1211-I-1437 is the first 
unbound subdwarf B star, for which we can rule out the slingshot mechanism. Alternatively, we may assume that the 
stars are old population II stars and therefore have to be bound. Then the kinematics of J1211+1437 set a lower limit 
of 2 X 10^^ M0 to the mass of the Galactic dark matter halo. 

Key words, stars: kinematics - stars: subdwarfs - stars: atmospheres - Proper motions - Galaxy: halo - line: profiles 



J> 1. Introduction Globular clusters, satellite galaxies, as well as large sam- 
r pies of halo stars, respectively, have been used to estimate 

^ The properties of the dark matter halo are important to un- j^^j^ mass. Actually only the objects with the most ex- 

H derstand how the Galaxy formed and evolved. Observations ^^^^^ velocities provide tight constraints and, hence, the 

. . ^ . of halo stars put constramts on theoretical models of halo j^^ss estimates depend mostly on them ([Sakamoto et all 

formation (e.g. [ Navarro et al. || 1996D . Lar^e surveys, such as 200I ISmith et al.llSh . A halo mass of about 2 x 10^^ 

the Sloan Digital Sky Survey (SDSS, [York et al.||2000|) and favoure d in earlie r investigations (iWilkinson fc Evans! 
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mm , provide large numbers of stars to trace the halo prop- f^^. b wer masses o f abo ut half tha t value (iBattadia et all 

erties, such as the total mass of the halo. |^0| ISmith et al.ll2007t iXue et al.l[200l . 



The hyper-velocity stars (HVS iBrown et al] 120051: 
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Kinematical studies of the hyper-velocit y stars were 
based on their radial velocities only. Recently. iTillich et alJ 
()2009f) were able to measure proper motions of an A-type 
HVS and study its 3-D kinematics to trace its place of 
birth in the Galactic disc. They found it to originate far 
from the Galactic centr e, thereby ch allengin g the SMBH- 
sling shot mechanism of iHillsl (jl988[ ) . Hence ITillich et aD 
(|200 9) suggested a runaway mechanism for the star's forma- 
tion. Further evidence that such a mechanism works comes 
from tw o similar studies o f the hyper run-a way stars HP 
217791 ()Heber et al.ll2008D and HIP 60350 ()Irrgang et al.l 
[20T0[ ). which were also found to originate in the outer rim 
of the Galactic disc nowhere near the Galactic centre. 



While most of the 17 HVS known today (jBrown et al.l 
I2009t ITillich et all I2009D are early-type ' mam-sequence 
stars, there is just one evolved low-n iass star, US 708, a hot 
subdwarf star of spectral type sdO (|Hirsch et al.ll2005l ). 

Most of the previous studies of halo stars to constrain 
the dark matter properties are hampered by the lack of 
proper motion measurements. Hence they had to rely sub- 
stantially on radial velocity distributions. In such cases only 
four coordinates (i.e. two position values, distance and ra- 
dial velocity, RV) of the 6D phase space are determined 
and the missing proper motion components are handled in 
a statistical appr oach. In the presently most extensive study 
IXue et al] (12008') measured radial velocities for more than 
10,000 blue halo stars from the SDSS and classified their 
sample as a mix of blue horizontal branch (BHB) stars, blue 
stragglers and main-sequence stars with effective tempera- 
tures ro ughly betwe e n 7,00 and 10,000 K according to their 
colours. IXue et al.l (j2008[ ) selected 2400 blue horizontal- 
branch stars to estimate the halo m ass out to 60 kp c to be 
1.0 X 10^^ Mq using a halo mo del of Navarro et al.l 



pc to oe 
111991 . 



For one star from that sample iPrzvbilla et al.l ( 201C ) were 
able to obtain proper motion and carry out a detailed kine- 
matic analysis, which revealed an inbound Population II 
horizontal branch star with a Galactic rest-frame (GRF) 
velocity of ^700 km s^^ at its current position. This makes 
it the fastest halo star known, and provided a lower limit 
of 1.7 X 10"'^^ M0 for the total halo mass of t he Galaxy, 
sign ificantly exceeding the value determined by IXue et al] 



sign il 
(|200? 



2008). 

This example shows that it is rewarding to study the 
kinematics of additional stars in the halo and to con- 
sider classes of stars other than BHB stars, as well. Of 
course, the Galactic h alo hosts a plethora of white dwarfs 
(|Eisenstein et al. ' 2006). However, they are so faint that 
they can be analysed in the solar neighbourhood only. 
Another group of evolved low mass stars are the hot subd- 
warf stars (sdB, sdO) that dominate the population of faint 
blue stars at hig h Galactic latitu des to visual magnitudes of 
about V=18 (Gr een et al.lll986l ). They are considered to be 
helium core burning stars with very thin (<O.O2M0) inert 
hydrogen envelop es and masse s around 0.5 Mq. Following 
ideas outlined bv iHebeJ (|l986l ). the sdBs can be identified 
with models for extreme horizontal branch (EHB) stars. 
An EHB star bears great resemblance to a helium main- 
sequence star of half a solar mass and it should evolve sim- 
ilarly, i.e. directly to the white dwarf cooling sequence, b; 
passing a second giant phase (for a review see iHeber 2009 
For the for mation of subdwar f B stars three scenarios arc 
discussed bv lHan et ahl (|2003D : common envelope ejection, 
stable Roche lobe overflow (RLOF) , and the merger of two 



helium white dwarfs. Some alternate sce narios for the for - 
mation of single sdB stars are reviewed bv l0stensenl (|2009l ). 

Hot subdwarf stars exist in the field of the Galaxy 
but also in globular clusters, in the Galactic bulge and 
have even been re solved in the elliptical galaxy M 32 
(iBrown et al.l 20081). K inematical studies (|Altmann et al.l 
12004 ; NaDiwotzkill2008f ) indicate that they occur in all stel- 
lar populations of the Galaxy. 

However, very little is known about the halo population 
of hot subdwarfs except those in globular clusters. Some 
high- velocity hot subdwarfs have attracted interest because 
of their high radial velocities, most notably, the sdO star 
US 708, whose radia l velocity in the r est- frame was mea- 
sured at 751 km s~^ (|Hirsch et al.ll2005D - the second HVS 
star discovered. Unfortunately we cannot deduce the origin 
of the star, as we lack a reliable proper motion measure- 
ment. 

Motivated by the discovery of US 708, we embarked 
on a project to identify a sample of population II hot 
subdwarfs and study their kinematics from radial veloc- 
ity and proper motion. W e make use of the MUCHFUSS 
survey (|Geier et al.ll20l6b[l . which searches for close bina- 
ries with high radial-velocity variations. The search strat- 
egy also provides targets that are not close binaries but 
travel through space at high RV without variations. These 
stars are the targets of our investigation. Accordingly we 
entitled our project Hyper-MUCHFUSS as we provide an 
extension of MUCHFUSS. 

The paper is organised as follows. In Sect. |2|we intro- 
duce our survey for HVS and in Sect. |3| our sophisticated 
proper motion measurement method. The kinematical ana- 
lysis techniques are shown in Sect. Hi In Sect. [5] we present 
our results and summarise and conclude in Sect. [5] 



2. Survey 

The enormous SDSS database is well evaluated in terms 
of errors and accuracy. Hence it is the perfect start- 
ing point for the MUCHFUSS survey. In order to se- 
lect subdwarf candidates, we used several indicators, e.g. 
colour, spectral classification and radial velocity (RV). 
Fig. [T] shows a fiowchart of the target selection method. 
We select ed sdO/B c andidates by colour (g-r<0.1 and u- 
g<0.4, see iGeier et al, 2010b) and pre-classified their spec- 
tra by visual inspection. Measuring the radial velocity 
by fitting synthetic models, we selected only stars faster 
than ±100 km s~^. For most of the known s dB binaries, 
the RV semi- amplitudes are below 100 km s~^ (|Geier et al.l 
l2010bt ). Hence the radial velocity of a typical sdB binary 
of the Galactic disc will rarely exceed 100 km s^^ in ab- 
solute value and these stars are consequently excluded in 
MUCHFUSS. We converted the heliocentric RV to the 
Galactic rest-frame (GRF). The larger the RV the higher 
is the priority we assign to the target for our survey. 
Especially stars with absolute GRF velocities of more than 
275 km s~^ are high priority HVS candidate s. Thi s obser- 
vational cut was introduced by Brown et al] (|2007l ) to dis- 
tinguish HVS from halo stars by their origin. More infor- 
mation on MU CHFUSS and the pr ocess of target selection 
is presented bv lGeier et al.l (j2010b[ ). 

More than 250 targets with measured radial velocity 
remained on the MUCHFUSS target list, serving as a 
first epoch. Second-epoch spectra were obtained with ESO- 
VLT/FORS (i? w 1800, A w 3730 - 5200A), WHT/ISIS 
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Fig. 1. Target selection method based on SDSS. Stars are 
selected for which the absolute value of the RV exceeds 
100 km s~^. Hence disc stars are efficiently removed. 



{R w 1800, A w 3730 - 5200A), CAHA-3.5m/TWIN {R w 
4000, A « 3460 - 5630A) and ES0-NTT/EF0SC2 (i? w 
2200, A « 4450 - 5110A). Due to the different wavelength 
coverage and resolution, the number of visible absorption 
lines is correspondingly restrictecfl. Up to now, such spectra 
are available for 88 stars, which is about one third of the 
target list. The RVs have been measured by x^-fitting of 
suitable synthetic spectra. Based on the MUCHFUSS RV 
list, we regard a star as RV constant if its velocity is consis- 
tent with the first epoch within the respective error limits. 
However it is obvious that these errors individually depend 
on a variety of conditions like e.g. the S/N (see Table[T]), the 
resolution, wavelength coverage and the number of visible 
absorption lines. 

The MUCHFUSS results for the b inaries will be pre- 
sented in a dedicated paper bv lGeier et al. (2010a, m prep.j. 
While 49 of the target stars were found to be RV variable, 
39 do not show variations and therefore make up the sam- 
ple that is studied in this paper. In the following sections 
we will exclusively focus on these candidates. 



3. Proper motion 

A kinematical analysis can only be done, if the star's lo- 
cation and space velocity is known. Hence we attempt to 
measure proper motions and determine spectroscopic dis- 
tances for all our 39 candidates. We collected all avail- 
able independent position measurements on Schmidt plates 



^ Note that the EFOSC2 spectra cover only H^j and 2 He I 
lines. 
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Fig. 2. Proper motion components derived from the posi- 
tion measurements for J1211+1437, where 1978.77 is the 
zero epoch. 



(APM - iMcMahon et"all[2000t SSS - iHamblv et ahllMl 
and combined t hem with the SDSS and other avai lable po- 
s itions (CMC14 ICarlsberg-Meridian-Ca talog''2006': 2MASS 
- iCutri et al.ll2003t UKIDSS - iLawrence et al. 2007) . Multi- 



epoch SDSS astrometric measurements, taken from the 
Princeton SDSS database (jFinkbeiner et al.ll2004[ ) were in- 
cluded in the proper motion fit. We obtained more mea- 
surements of Schmidt plates, from up to 14 different epochs 
in case of overlapping plates of the Digitised Sky SurvejQ 
(DSS). The FITS images of 15 by 15 arcmin size were ex- 
tracted from all available plates and the ESO MIDAS tool 
center/gauss was used to measure positions. To measure ab- 
solute proper motions, we initially had to find and identify 
compact background galaxies from the SDSS. It is highly 
important to obtain enough galaxies well distributed over 
the image, although galaxies tend to cluster. We used at 
least 10 galaxies per measurement and determined the posi- 
tion of the star relative to the reference galaxies. Finally we 
transformed the target positions on all the Schmidt plates 
to the SDSS system. The small fields allowed us to apply 
a simple model (shift-f rotation) . The typical error in both 
coordinates at one epoch, including the error of the trans- 



http:/ /archive. stsci.edu/cgi-bin/dss_plateJinder 
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Table 1. Heliocentric radial velocity table of the target sample with SNR per pixel in the continuum next to H^. 



name 






short name 


type 


OBS 


MJD 




Vrad 

1 —1 

km s 


bNR(Q)H^ 


SDSSJ084447, 


,00+113910. 


,0 


J0844+1139 


sdB 


SDSS 


2453801. 


.1251331 


+202.7 ±5.0 


52 












NTT 


2454755. 


.86312183 


+219.5 ±30.4 


13 


SDSSJ084556, 


.85+135211. 


.3 


J0845+1352 


sdB 


SDSS 


2453799. 


.765353 


+115.1 ±9.6 


33 












NTT 


2454756. 


.86725607 


+ 92.0 ±15.1 


22 


SDSSJ094850, 


.47+551631. 


.7 


J0948+5516 


sdB 


SDSS 


2451991. 


.6335069 


-138.3 ±12.2 


17 












WHT 


2454587. 


.46586414 


-158.5 ±17.8 


9 












WHT 


2454588. 


.38885923 


-127.4 ±6.6 


15 


SDSSJ102057, 


.16+013751. 


.3 


J1020+0137 


sdB 


SDSS 


2451997. 


.186412 


+243.3 ±8.3 


48 












FORS 


2454638. 


.47970756 


+239.4 ±7.8 


36 


SDSSJ121150, 


.27+143716. 


.2 


J1211+1437 


sdB 


SDSS 


2453466. 


.7114641 


+220.2 ±9.5 


32 












FORS 


2454575. 


.76860445 


+207.5 ±25.6 


11 


SDSSJ155635, 


.80+470851. 


.8 


J1556+4708 


sdB 


SDSS 


2452354. 


.493235 


-393.5 ±19.9 


18 












WHT 


2454586. 


.60330613 


-385.4 ±7.4 


14 


SDSSJ163213, 


.05+205124. 


.0 


J1632+2051 


sdB 


SDSS 


2453224. 


.6422685 


-238.6 ±6.7 


38 












FORS 


2454656. 


.62186841 


-237.9 ±2.8 


45 


SDSSJ164419, 


.44+452326. 


.7 


J1644+4523 


sdB 


SDSS 


2452051. 


.9386574 


-314.1 ±4.5 


32 












WHT 


2454586. 


.71567501 


-312.8 ±11.6 


31 


SDSSJ215648, 


.71+003620. 


.7 


J2156+0036 


sdB 


SDSS 


2452933. 


.5838194 


-177.2 ±6.5 


32 












NTT 


2454757. 


.70342858 


-190.1 ±37.2 


13 


SDSSJ224451, 


.81+010630. 


.9 


J2244+0106 


sdB 


SDSS 


2452146. 


.8122801 


-166.0 ±5.9 


24 












NTT 


2454756. 


.61349441 


-149.1 ±19.2 


17 


SDSSJ084938, 


.85+145503. 


.2 


J0849+1455 


DA 


SDSS 


2453818. 


.6548958 


+119.0 ±37.2 


15 












NTT 


2454757. 


.87125523 


+ 82.6 ±33.9 


9 


SDSSJ135824, 


.61+472931. 


.6 


J1358+4729 


DA 


SDSS 


2452723. 


.9695428 


-260.6 ±14.2 


24 












WHT 


2454586. 


.54874204 


-257.0 ±14.3 


22 



Table 2. PM measurements for the positive detections to- 
gether with the number of epochs riep and the number of 
reference galaxies Ugai- The V magnitudes have been de- 
rived from SPSS photometry using the transformations of 
iJordi et alJ \200Q . 



name 


V 








{^a cos (5)) 




Hep / ^gal 




mag 






mas/yr 


mas/yr 




J0844+1139 


16. 


.11 


±0, 


.02 


-11.5 ±2.6 


-11.9 ±2.1 


6/ 7 


J0845+1352 


17. 


.42 


±0, 


.02 


-10.7 ±2.3 


-3.7 ±2.6 


5/14 


J0948+5516 


18. 


.59 


±0, 


.03 


-1.9 ±2.9 


-7.4 ±1.7 


7/12 


J1020+0137 


16. 


.91 


±0, 


.02 


-2.4 ±1.3 


-7.3 ±1.2 


12/15 


J1211+1437 


17. 


.87 


±0, 


.03 


-12.1 ±1.8 


-27.2 ±1.4 


13/14 


J1556+4708 


18. 


.61 


±0, 


.03 


-6.3 ±1.9 


+1.3 ±1.8 


10/14 


J1632+2051 


17. 


.62 


±0, 


.03 


-12.5 ±3.0 


-1.6 ±3.6 


7/16 


J1644+4523 


17. 


.36 


±0, 


.02 


+4.7 ±2.8 


-26.1 ±3.3 


8/13 


J2156+0036 


17. 


.91 


±0, 


.03 


-1.3 ±1.6 


-7.4 ±1.2 


26/13 


J2244+0106 


18. 


.14 


±0, 


.02 


+2.1 ±0.8 


-3.7 ±1.3 


30/14 


J0849+1455 


18. 


.85 


±0, 


.03 


-49.5 ±1.8 


+6.4 ±1.6 


8/ 9 


J1358+4729 


18. 


.03 


±0, 


.03 


+76.1 ±2.4 


+ 14.4 ±2.4 


14/15 



formation using the measured galaxies and the individual 
centroiding error for the target, was 150mas. We obtained 
one position per epoch and used linear regression to derive 
the proper motions with their errors (see e.g. Fig [2]). 

Especially for the old POSS-I epoch we noticed large 
discrepancies between the positions, which are probably 
due to the different colour filters and quality of the plates. 
The best way to minimize possible systematic effects is sim- 
ply to use the same set of reference galaxies for every epoch. 



We regard a proper motion as detected, if the position mea- 
surements do not show a large spread relative to the linear 
fit and the derived proper motion is significantly different 
from zero. 27 stars turned out to have a proper motion con- 
sistent with zero, while for 12 of our 39 candidates (w 38%) 
a significant proper motion was measured. One of them is 
the sdB J2156+0036 which has already been used to show 
the po tential of our method (see Fig. 1 of iTillich et al.l 
l2010a[ ). This star represents also a typical star with proper 
motion from our target sample. Note that for all but two ob- 
jects (J1644+4523 and J0849+1455) the determined proper 
motions agree within their er rors with those of t he recently 
published PPMXL catalogue (jRoeser et al.l20l"o[l , where for 
all objects our proper motion errors are smaller than the 
PPMXL ones. 

Tables[T]and[2]list the measurements for the sample ana- 
lysed in this paper. In Fig [5] we demonstrate the measure- 
ment of the proper motion for J1211+1437, an outstanding 
target that will be discussed in detail in section [Ql 

4. Distance and Kinematics 

The second important parameter is the distance, necessary 
to determine the current location of the stars and their 
transversal velocities along with their proper motions. We 
determine the distance from a quantitative spectral anal- 
yse, which provides the effective temperature, gravity and 
helium abundance. The distance is derived from the atmo- 
spheric parameters and the apparent magnitude, by adopt- 
ing the canonical mass of 0.48 Mq. 

In most of the cases the stellar spectra from the SDSS 
provide sufficient S/N and resolution for a quantitative 
spectral analysis. Furthermore, in some cases they make up 



4 



A. Tillich et al.: The Hyper-MUCHFUSS project: probing the Galactic halo with sdB stars 



the most reliable data we have. We applied x^-fitting of syn- 
thetic line profiles to the Balmer and helium lines in order 
to determine the atmospheric parameters and abundances 
(see Fig. Inevitably, our selection procedure produces 
some misclassifications, therefore we also ended up with two 
DA white dwarfs in our sample. All of the subdwarfs have 
been analysed using fully metal line-blanketed LTE mod- 
els of solar metallicity (H eber et al. 2000j3- Finally spec- 
troscopic distai ices are calculated usin g the astrophysical 
fluxes following iRamspeck et al.l (|200l[ ). Independent stud- 
ies showed, that the atmospheric parameters (and hence the 
distances) depend only little on the choice of the metallicity. 
iHeber et al.l (|2000[ ) studied the effect by comparing results 
from models with solar composition to those of 1/100 solar 
metallicity in a detailed study of high resolution spectra 
for three sdB stars. The differences are of the order of 200- 
300 K in Toff and 0.03 dex in \ogg, far lower than the 
uncertainties of our results. Hence, metallicity effects are 
not significant. The uncertainties of the distance have been 
determined usi ng a Monte Carlo a.pproach for the input 
quantities (cf. iTiUich et al.) l2009l l2010bl: iPrzvbilla et al] 
120101) . 

Ap plying the Galactic potential of lAllen fc Santillaiil 
()1991|) we calculated orbits and reconstructed the path 
of the star with the program of lOdenkirchen fc Broschd 
(|1992D . The distance of the GC from the Sun was adopted 
to be 8.0 kpc and the Sun's motion w ith respect to the 
local standard of rest was taken from iDehnen fc BinnevI 
(119981) . Based on our data we are not able to draw con- 
clusions on detailed structural properties of different galac- 
tic potentials. However, we are able to discuss global ef- 
fects (see Section [5]) . The error in space motion is domi- 
nated by that of the distance (via the gravity error) and 
those of the proper motion component^. Varying these 
three quantities within their respective errors we applied 
a Monte Carlo procedure for the error propagation and to 
derive the median GRF velocities at t he present location 
and their distribution (see Fig. El cf. iTillich et al1l2009l 
l2010bHPrzvbilla et al.|[20Tol ). We compared our results with 
the local esc ape velocities as c a lculat ed from the Galactic 
potential of lAllen fc SantillanI (|1991[ ). Furthermore, from 
kinematic characteristics (U, V, eccentricity e and the z- 
component of the angular momentum 3z) we obtained ad- 
ditional information about the population membership of 
the stars. Following the trajectories backwards in time we 
can even put constraints on the origin of the stars, applying 
the same Monte Carlo method. 

5. Results 

Until now about one third of the targets of the MUCHFUSS 
survey have been observed. Accordingly 39 of them have 
constant RV within the detection limits. Twelve stars 
showed proper motions significantly different from zero 
and have been analysed in detail using quantitative spec- 
troscopy. In Fig. |6] we show a comparison between the ob- 
served and the synthetic spectra for two subdwarf stars. 



For t he DA white dwarfs we used the synthetic models by 
iKoested (2009). 

* Although some of the second-epoch spectra are of low qual- 
ity, the RV errors (see Table [T]) are irrelevant for the error bud- 
get. The distribution of the space velocity is dominated by the 
errors in the transversal velocity components, which may have 
uncertainties as large as ±100 km s~^. 




Fig. 3. Comparison of the position of the 10 subdwarfs 
(with error bars) from HMF p roject in the Teff-lo g q- 
diagram to evolutionary tracks of iDorman et al.l (1 19931) for 
a metaUicity of [Fe/H]^ -1.48. The ZAHB and the TAHB 
are indicated only from fitting of the tr acks. The SPY sub - 
dwarf B sample is plotted for reference (jLisker et al.ll2005h . 
Note that two subdwarfs lie lie well above the canonical 
TAHB and therefore most likely are in a shell-He burning 
post-EHB stage of evolution. 



Table 3. Stellar parameters for the remaining stars with 
PM significantly different from zero. 



name 


type 


Tcff 
K 


logff 


og nHc/nn 


J0844-hll39 


sdB 


28600 


±740 


5.33 


±0.12 


-2.6 ±0.4 


J0845-H352 


sdB 


24570 


±1070 


5.59 


±0.17 


<-3.0 


J0948-)-5516 


sdB 


33890 


±560 


5.84 


±0.20 


-1.5 ±0.2 


J1020+0137 


sdB 


28960 


±760 


5.39 


±0.14 


<-3.0 


J1211+1437 


sdB 


32530 


±580 


5.52 


±0.16 


-2.9 ±0.5 


J1556+4708 


sdB 


34410 


±1280 


5.34 


±0.25 


<-3.0 


J1632-f2051 


sdB 


26870 


±610 


5.31 


±0.09 


-2.1 ±0.2 


J1644-h4523 


sdB 


31680 


±410 


5.78 


±0.11 


-2.9 ±0.3 


J2156+0036 


sdB 


29080 


±810 


5.60 


±0.15 


-2.5 ±0.3 


J2244+0106 


sdB 


33580 


±680 


4.75 


±0.20 


-1.5 ±0.2 


J0849-hl455 


DA 


27750 


±830 


7.32 


±0.15 




J1358-I-4729 


DA 


10640 


±110 


8.10 


±0.07 





Two stars turned out to be DA white dwarfs rather than 
sdB stars. The derived stellar parameters are shown in 
Table [31 while the respective RV can be found in Table HI 
In Fig. Owe show a Toff-log g-diagram for the 10 sub- 
dwarfs in comparison to the refere n ce sam ple of the SPY 
sdB stars analysed bv iLisker et al.l (|2005[ ). Obviously 8 of 
our subdwarf B stars reside on the extended horizontal- 
branch (EHB) within the respective errors. However, the 
stars J1556+4708 and J2244+0106 lie well above the EHB, 
and must therefore be in a post-EHB stage of evolution 
if they have the canonical mass of 0.48 M©. Fig. [3] shows 
that our subdwarf B sample goes well together with the 
SPY sample. Based on these parameters we determined 
the distances, which are given in Table HI Together with 
the median space velocity we als o obtained the escape ve - 
locity in the Galactic potential of lAllen fc SantillanI ()1991l ). 
To further quantify the kinematics of the stars, we made a 
comparison with the kinematics of white dwarfs. A sample 
of 398 white dwar fs from the SPY survey was studied by 
iPauli et all (|2006[ ). Based on the 3D-orbit, the V - U di- 
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Table 4. Estimated distance, escape and GRF velocities 
(median) and loc al escape velocities for the Galactic po- 
tential of I Allen fc Santillan (1991). 



name 


dist 




VGRF 






kpc 


km s 


km s 


km s 


J0844+1139 


3.14 ±0.26 


202.7 ±5.0 


2021^^ 


513 


J0845+1352 


3.50 ±0.38 


115.1 ±9.6 


230t^° 


511 


J0948+5516 


6.51 ±0.88 


-142.1 ±8.6 


106t^^ 


490 


J1020+0137 


4.31 ±0.39 


238.0 ±8.3 


107+13 


513 


J1211+1437 


6.45 ±0.68 


220.1 ±9.5 




507 


J1556+4708 


11.76 ±2.03 


-385.4 ±7.4 


qt:q+174 


487 


J1632+2051 


6.09 ±0.36 


-236.1 ±3.2 


97r + 120 
^' J_110 


544 


J1644+4523 


3.68 ±0.27 


-314.1 ±4.5 


468l^r 


534 


J2156+0036 


5.37 ±0.54 


-177.2 ±6.5 


fic+58 
80_44 


534 


J2244+0106 


18.29 ±2.45 


-166.0 ±5.8 




458 


J0849+1455 


1.11 ±0.11 


119.0 ±37.2 


354+g 


531 


J1358+4729 


0.14 ±0.01 


-260.6 ±14.2 


295+1? 


541 



agram and the e- Jz-diagram they introduced a kinematic 
population classification scheme and combined it with age 
information. They performed a detailed kinematical analy- 
sis accounting for errors by means of a Monte Carlo error 
propagation code, similar to our method. A sample of abun- 
dance selected typical main-sequence stars served as refer- 
ence sample. They derived the Scr-contours for the V — U 
diagram on which the kinematic classification is based. For 
the e- Jz diagram a "Region B" is defined such that it ex- 
cludes as many thick-disc stars as possible. The last crite- 
rion is the p-Z-diagram, which is used to classify the or- 
bits among the populations by c omparing with template 
Galactic orbits ()Pauli et al.ll2003[) . A substantial thick disc 
fraction of 7% was found, while only 2% of the DAs show 
characteristic halo properties. 

In Fig , m we compare ou r sdB sample to the white dwarf 
sample of lPauh et all()2003D . All sdB stars lie far away from 
the thin disc population of the white dwarfs. They are found 
in those regions of both diagrams where very few white 
dwarfs lie; only the rare ones belonging to the halo popula- 
tion. The sdB J0845+1352 shows thick disc kinematics (see 
Fig. |4]). Nine sdB stars possess halo characteristics, as they 
reside clearly outside the Str thick disc contour in the dia- 
grams. According to the e- Jz diagram, the sdBs can be di- 
vided into two subgroups. Furthermore, the 3D-orbits (see 
Fig. [S] for examples of an Gl and a G2 star, respectively) 
are all clearly favour a halo membership. However, for 10 
stars at least two of three indicatio ns are present, w hich 
renders them as halo stars, following iPauli et al.l ()2006[ ). In 
the e- Jz diagram the cluster of subdwarfs with high eccen- 
tricities at Jz w clearly catches one's eye. We call this 
cluster "group 1" (Gl), as the respective stars share simi- 
lar kinematic properties. They have nothing to do with the 
disc rotation and show only marginal velocity components 
in the direction of the Galactic plane. Such a behaviour is 
very typical of halo objects, as they travel mainly perpen- 
dicular through the disc. Four subdwarfs have very negative 
angular momenta, which means that they are on retrograde 
orbits. This group of subdwarfs we term "group 2" (G2). 



The two targets with the highest measured PMs turned 
out to be DA white dwarfs. In fact, J1358+4729 was discov- 
ered and analysed already by Eisenstein et al. (2006). They 
derived atmospheric parameters (Tcff — 10635 ±61 K, logg 
— 8. 16 ±0.053), which are perfectly consistent with our val- 
ues. J1358+4729 shows thick disc kinematics (see Fig. S]) 
according to both diagrams and does not belong to the 
halo population. The other white dwarf J0849+1455 has 
no t been analyse d before, except for a proper motion study 
by lEvand (|l992f) , and belongs to the halo population (see 
Fig. U]). In their last orbits their trajectories never came 
close to the central part of the Galaxy. 

5.1. J1211+1437: extreme halo or hyper-velocity star 

According to the kinematic analysis presented above the 
sdB star J1211+1437 shows the most extreme kinematics 
in G2 of retrograde orbit stars, as it lies far away from the 
other stars in the kinematic diagrams. 

J1211+1437 - a HVS candidate 

The Galactic rest-frame velocity of J1211+1437 is so high 
(vGflF = 713^139 km s~^, see Fig |5]) that it exceeds the 
local Galactic escape velocity o f Veae = 507 km s~^ if w e 
adopt the Galactic potential of lAUen fc SantillanI (|1991[ ). 
J1211+1437 is a HVS candidate that could have been 
ejected from the Galactic centre by the Hills mechanism. 
In order to test this hypothesis we traced the trajectory 
back to zero Galactic latitude. 

As can be seen from Fig. [7] J1211+1437 does not orig- 
inate in the GG. Its place of origin is more likely to lie in 
the Galactic dis c, where no SMBH is known to exist. Hence, 
the iHilld ()1988l) slingshot mechanism can by excluded. 

The ejection velocity for that star, if it would have be- 
longed to a co-rotating disc population, is about Vej ~ 
550 km s~^. Hence the star might have been ejected from 
the Galactic disc as a run-away star. 

J1211+1437 and the dark matter halo 

If, however, J1211+1437 is an old halo star that experi- 
enced several disc passages the star has to be bound to 
the Galaxy. This can only be achieved if the Dark Matter 
halo is more massive than the adopted one in the poten- 
tial of lAllen fc SantillanI (jl99lt) . In general the kinemat- 
ics of the most extreme stars p rovide the best constraint s 
on the mass of the d ark halo ([Wilkinson fc Evans! [19991 ) . 
iPrzvbilla et al.l (j2010l) found a population II BHB star trav- 
elling towards us at such a high speed that a dark mat- 
ter halo mass of at least 1.7 x 10^^ Mq is required to 
keep the star bound to the Gal axy. This value is consid- 
eral3ly larger than adopted by lAllen fc SantillanI ()1991[ ). 
Therefore numerical experiments were carried out in which 
the Galactic potential was modified by increasing the mass 
of the dark matter halo. This constrains the total Galactic 
mass to exceed M^^^^i = 2.0±tli x 10^^ Mq. The errors are 
based on the Monte Garlo distribution of the space veloc- 
ity. Although the uncertainties are large, our value is per- 
fectly consisten t with other pro per motion based kinematic 
mass estimates (iPrzybilla et al.. 2010: Sakamoto ctjal,200i 
IWilkinson fc Evanslll999ir which also support a high-mass 
Galactic halo. However, it has to be stressed that this con- 
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Fig. 4. V-U (left) and e-3z diagram (right) for our 12 targets. The white dwarf sample (+) of iPauli et alJ (|2006[ ) serves 
as reference. The solid ellipses render the Scr-thin and thick disc contours, while the solid box marks the thick disc region 
(Region B). Note that 10 of our targets are outside of the regions for the two disc populations in both diagrams, but 
form two distinct groups. The 5 stars inside the dashed ellipse are called group 1 (Gl), while the stars with negative 3z 
are defined as group 2 (G2). 





p / kpc p / kpc 

Fig. 5. p-Z-diagrams for the two group representatives, J1020+0137 (group 1, left) and J1632+2051 (group 2, right). 



elusion can only be drawn, if the orbit of J1211+1437 is 
bounc0. 

5.2. J 1644+4523- the fastest Gl object 

Another star that is worth closer inspection is the sdB 
star J1644+4523, which belongs kinematically to the Gl 
group. The star is heading towards us and therefore must 
be bound. It is the second fastest hot subdwarf of the 
halo population. We measured a GRF velocity of vgrf = 
468^gg^ km s~^ , which is below the escape velocity of Vesc = 
534 km s^^ in the Galactic potential of I Allen fc SantillanI 
(|l99lh . As for J1211+2423 we might consider alternatively 
that the star was ejected as a run-away star. We calcu- 
lated the trajectory backwards into the past and found that 

^ The rest of the sample is bound to the Galaxy, al- 
though the subdwarfs J1556+4708 and J1644+4523 might ex- 
ceed the Galactic escape velocity if their respective errors (see 
Table [H) are taken into account and the Galactic potential of 
I Allen fc SantillanI (|199ir ) is applied. 



the star might have originated in the Galactic bulge (GB), 
which can be seen from Fig. [HI The last pericenter pas- 
sage occurred at a distance of only ~ 0.5 ± 1.0 kpc from 
the Galactic centre (see Fig. [7]) and the apocenter distance 
of the star's trajectory is located far out in the halo at 
^ 115 kpc. The high ejection velocity of Vej = 597 km s^^ 
could possibly be consistent with the ejection by the SMBH 
in the GC, as the error ellipse encloses the GC (see Fig. [7]). 
The time-of-flight (TOF) is much longer [TOF = 1.27 Gyr) 
than the helium core burning phase lasts (w 200 Myr). This 
means that J1644-I-4523 evolved to a subdwarf a long time 
after it was possibly ejected from the GB. This is very prob- 
lematic, as EHB stars most likely requir e a binary interac- 
tion scenario to form ()Han et al.l |2003[ ). In this case the 
ejection by the SMBH would disrupt any initial binary and 
the sdB cannot form. One solution could be a shorter TOF. 
The kinematic analysis of J1211-t-2423 indicated that the 
dark matter halo may be more m assive th an ass umed in the 
Galactic potential of lAUen fc SantillanI (|l99lh . Therefore 
we repeated the kinematical calculations for a higher dark 
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Fig. 6. Comparison between the synthetic model with the observed spectrum (WHT/SDSS) for the sdBs J1644+4523 
(top) and J1211+1437 (bottom). The overall agreement is good. 



halo mas s and adopted Mi^! Z = 3.4 x lO^^ M© as sug- 
gested bv lAbadi etall (I2009D with a mass distribution out 
to 100 kpc. following Ullen fc SantillanI ()199H ). The place 
of the last disc intersection remained almost unchanged but 
the time-of-flight is strongly reduced to TOF = 145 Myr, 
which is consistent with the lifetime of hot subdwarfs. 
Hence, in a high-mass halo the star could have been ejected 
as subdwarf and did not have to evolve after the ejection. In 
that case the last pericenter passage occurred at a slightly 
higher distance of ~ 0.65 kpc from the Galactic centre and 
the apocenter distance of the stars trajectory is located in 
the halo at only ^ 22 kpc. 



5.3. The entire sdB sample 



Hence we consider it worthwhile to inspect the run-away 
scenario for all target sdB stars by taking a closer look at 
their trajectories from the Galactic plane to their present 
location in the Galactic halo. Fig. [TT] displays the areas in 
the plane from which the stars could have originated for the 
kinematic group Gl (top) and G2 (bottom). As can be seen 
the Gl stars would originate from the inner Galactic disc 
or bulge, in particular J1644-t-4523, whereas the G2 stars 
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Fig. 9. 3D plot of the current positions of our 12 stars (black triangles) relative to the Galactic Disk (grey), showing also 
the trajectory for J1644+4523 and the area of potential origin (dark grey). 
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Fig. 7. Intersection region of the subdwarfs' past trajecto- 
ries for J1211-I-1437 and J1644-f4523 with the GD in order 
to conclude on their origin (grey region). For orientation 
the position of the sun is marked as well. 
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Fig. 8. Galactic restframe velocity distribution for 

J1211+1437 derived from a Monte Car lo pr ocedure 

wi th a depth of 10 000 (cp. iTillich et al.l [20091 l2010bl: 
iPrzvbilla et al.l [20101 ). The gray shaded area indicate the 
velocities for whi c h the star is bound in the potential of 
lAllen fc SantHIm] ([Tool . 



(except J1632-j-2051) would come from the outer dis<3. The 
case of J1211-t-1437 demonstrates that some of our targets 
could possibly be ejected run-away stars. If so the places of 
plane crossings of the Gl stars would not rule out an origin 
in the Galactic centre (except for J0120+0137). However, 



The remaining thick disc subdwarf J0845-I-1352 shows a tra- 
jectory, which intersects within the Galactic disc. 



a Galactic centre origin of the G2 stars is excluded as their 
potential places of origin do not include the central part of 
the Galaxy. 
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Fig. 10. Kinematic V-U diagram for the programme 
sdB stars in compa rison with the sample (+ signs) of 
Hit mann et all (|2Q04f ). The solid cUipscs render the Scr-thin 
and thick disc contours of Pauh et al- QOOg). 



5.4. Kinematics of hot subdwarfs 

Up to now, the only 3D kin ematical study of suc h sdBs 
in the field was performed by lAltmann "eTaLl (l200l . From 
a sample of 114 sdBs they found 15 to have orbits, which 
differ considerably from disc orbits, while the vast major- 
ity is consistent with the disc population. lAltmann et alj 
(j2004l) noticed that their sample may contain close binary 
stars that are radial velocity variable, because the fraction 
is known to be about 40%. At the time of writing several 
stars of their sample had not been checked for radial veloc- 
ity variability, but have been shown to be RV variable later 
on. We do not correct for this effect. 

We compare the kinematical properties in the V-U- 
plane of t heir sample t o our s in Fig. 1101 Applying the 



criteria of Pauh et al.l (|2006[ ) the vast majority of the 



lAltmann et al. (2004) sample belong to the thin disc, with 
about two dozen thick disc and a dozen halo stars. The 
most eye-catching difference between the two samples in 
Fig. [To] is the absence of stars with rapid retrograde orbits, 
i.e . our class G2. Only very few halo sdBs in the sample 
of lAltmann et all ()2004[ ) would be classified as Gl (show- 
ing essentially no Galactic rotation). However, our sample 
lacks stars with large positive V velocities. 

When comparing the two sa mples we have t o con sider 
that the stars in the sample of lAltmann et al.l (I2004ri are 
much brighter and, hence, nearer by than those in our sam- 
ple. Consequently, we probe a much larger volume of space. 



6. Summary and conclusion 

Thirty-nine high-velocity faint blue stars have been iden- 
tified from the MUCHFUSS survey. We were able to mea- 
sure proper motions significantly different from zero for 12 
of them by determining positions of 60 years time base. 
Quantitative spectral analyses were performed to derive 
their distances. From an investigation of their kinematics 
nine sdB stars were identified as halo members and one 
as belonging to the thick disc. The two remaining stars 
turned out to be white dwarfs in the halo and the thick 
disc, respectively. Two distinctive kinematic groups among 



the sdB stars emerged: the normal halo subdwarfs (Gl) 
with low Galactic rotation and the extreme halo subdwarfs 
(G2) on highly-eccentric retrograde orbits. 

The extreme halo star J1211-I-1437 is of particular in- 
terest as it would be unbound to the Galaxy if the standard 
Galactic potential was applied. Hence it would be the sec- 
ond HVS among the hot subdwarfs and the first one of 
spectral type sdB. However, assuming the star is a popula- 
tion II object bound to the Galaxy, we can derive a lower 
limit to the mass of the dark matter halo of 2 x IO-'^^Mq. 
Other stars, J1644-I-4523 in particular, that are approach- 
ing the Earth have to be bound to the Galaxy and therefore 
might provide constraints on the dark matter mass. 

The existence of two kinematically distinct groups in- 
dicates different origins. In the ejection scenario the class 
Gl stars would originate from the inner Galaxy, whereas 
the G2 ones would come from its outer parts. However, if 
some of the sdB stars in our sample were run-away stars, 
we have to find an ejection mechanism. The Gl stars (ex- 
cept J0120-I-0137) may originate from the Galactic centre 
and the SMBH slingshot mechanism might work for them. 
For the G2 stars an origin in the Galactic centre can be 
excluded. Hence we have to find another mechanism to ex- 
plain their origin. 

Kinematic studies among massive B-type stars indi- 
cate that variation s of the typical runaway-scenario migh t 
provide an answer (|Heber et al.ll2008l: llrrgang et al.ll20100 . 
Gonveniently, sdBs mo st likely require a binary interac- 
tion scenario to form ()Han et al.l l2003f) . In principle one 
can think of a variation of th e RLO F mechanism as sug- 
gested bv iNelemans fc TaurisI (jlQQST ) in the context of the 
formation of undermassive white dwarfs. If the subdwarfs 
RGB progenitor is losing its envelope to a massive white 
dwarf, an asymmetric accretion induced collapse could oc- 
cur. The system can then be disrupted, with the accretor 
leaving t he system as a h i gh- ve locity neutron star. As com- 
puted bv lMarietta et al.l (|2000[ ). the companions SN explo- 
sion does not have to happen at the tip of the RGB for the 
envelope to be lost and the core to end up as an sdB star. 
In their calculations, the SN explosion itself can strip 96 to 
98% of the envelope from an RGB star, effectively leaving 
a naked He core. If the core is massive enough for helium 
burning it can experience a late core He flash, and end up 
as a single EHB star. As pointed out by Ostensen (20091), 
in both these cases the abandoned subdwarf would end up 
single and in an unusual gal ac tic or bit. 

Recently, IWang fc Hani (|2009D and lJustham et~aI1 
(2009) suggested a single degenerate SN la scenario to ex- 
plain the formation of the helium sdO US 708. Accordingly, 
a white dwarf accretes matter from am helium star in a 
close binary system. After exceeding the Chandrasekhar 
mass- limit, the white dwarf explodes as a type la super- 
nova. The orbital period has shrunk to about Ih and the 
corresponding orbital velocity of the hot He-sdO companion 
may exceed 500 km s^^. After the explosion of the primary, 
the binary is disrupted and the sdO companion is released 
at its orbital velocity. In addition it may be speeded up 
by the SN kick. However, all our programme stars are sdB 
stars, which retain a hydrogen-rich envelope of 10~^ to 10~^ 
Mq. Hence, it appears unlikely that the scenario applies to 
high-velocity sdB stars. 

As some globular clusters are known to host sdB stars, 
ejection from such clusters may give rise to high-velocity 
sdB stars. Alternatively, the programme stars may stem 
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from the disruption of a satellite galaxy in the Galactic 
halo. 

We shall extend our study as the MUCHFUSS survey 
proceeds. As demonstrated, the high-velocity hot subdwarfs 
are important tools to constrain the mass of the Galactic 
dark matter halo. Its full potential will develop once the 
GAIA mission will provide much more accurate distances 
and proper motions. 
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Fig. 11. Intersection region of the halo subdwarfs' past trajectories with the GD in order to conclude on their origin. 
Note that in the group Gl, 4 out of 5 subdwarfs are consistent with a GC origin (top). In the extreme group G2 we 
found 3 subdwarfs, which might originate in the outer Galactic rim (bottom) and have been likely produced by supernova 
runaway ejection. For orientation the position of the sun is marked as well. 
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